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Abstract 

After discussing the problem of lattice regularization of chiral gauge theories, a simple 
model for anomalous fermion number violation is formulated which can be numerically 
studied with present day technique. Exploratory results of numerical simulations of a 
■ two-dimensional U(l) Higgs model are presented. 

o 
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1 Introduction 

The anomalous baryon and lepton number violation is a consequence of the chirality of the 
electroweak gauge interactions M. This is a non-perturbative phenomenon and, therefore, a 
lattice formulation allowing for numerical simulations would be desirable. It is, however, well 
known that the non-perturbative formulation of chiral gauge theories is problematic (see, for 
instance, 0). 

In the first part of my talk some recent developments concerning lattice chiral fermions are 
reviewed and the possibility of existence of mirror fermions is discussed. Recent work on the 
limits of mirror fermion parameters is reviewed within the framework of a particular mixing 
scheme of fermion-mirror-fermion pairs. In the second part a simple model of anomalous 
fermion number violation is discussed which can be studied on the lattice ||. In the last part 
a simplified toy model in two dimensions, with chiral U(l)-symmetry, is discussed and some 
recent numerical results concerning topological charge and Chern-Simons number are presented. 



2 Chiral fermions on the lattice? 
2.1 Recent developments 

For an exact lattice formulation of the Standard Model on the lattice one has to solve the 
problem of chiral gauge theories. The recently proposed "domain wall fermion" formulation 
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[|], m can potentially lead to a solution by introducing an additional fifth dimension with four 
dimensional hypersurfaces representing the real world. A simpler, more economic, version of 
this idea is to throw away half of the five dimensional space, resulting in the "boundary fermion" 
formulation ||. The difficulty in these formulations arises at switching on the gauge coupling. 
If in the fifth direction the gauge coupling is chosen to be strong, one has to fight against the 
dangerous "layered phase" |7[], where subsequent four-dimensional hyperplanes do not interact 
and the theory becomes vector-like |§. In case of a four-dimensional gauge field, identical in all 
hyperplanes, the decoupling of the chiral fermion from its mirror partner with opposite chirality 
seems unlikely, at least as long as exact gauge invariance is maintained. 

A further obstacle for defining chiral gauge theories on the lattice is a recent extension 
of the Nielsen-Ninomiya no-go theorem including interactions and considering the continuum 
limit @. 

Avoiding the fermion doubling problem on random lattices has recently been also investi- 



gated, and does not seem to be promising, either JT^ . 

The remaining possibilities for the lattice formulation of chiral gauge theories are not ex- 
plicitly gauge invariant. In the modified Rome-approach capable to accomodate also fermion 



number violations [O] gauge fixing is essential. In case of the "reduced staggered fermions' 



12 1 gauge symmetry is explicitly broken by lattice actifacts, but hopefully restored by the 



dynamics. 

2.2 Mirror fermions 

Since the non-perturbative formulation of the electroweak interactions is so difficult, a natural 
question is whether "chirality" is perhaps only a low-energy phenomenon, and at high energy 
the space-reflection symmetry is restored by the existence of opposite chirality "mirror fer- 
mions" [jn|. If the presently known (almost complete) three fermion families were duplicated 
at the electroweak energy scale, in the range 100-1000 GeV, by three mirror fermion families 
with opposite chiralities and hence V+A couplings to the weak gauge vector bosons 0], then 



the whole fermion spectrum would be "vectorlike". This would very much simplify the non- 



perturbative lattice formulation of the Standard Model |L5|. 

The direct pair production of mirror fermions is not observed at LEP. This puts a lower 
limit on their masses of about 45 GeV. Heavier mirror fermions could be produced via their 
mixing to ordinary fermions. The present data imply some constraints on the mixing angles 
versus the masses which, however, strongly depend on the mixing scheme. (For an evaluation 
of the constraints implied by the high precision LEP data see the recent paper by Csaki and 
Csikor 111 
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2.3 Fermion— mirror- fermion mixing scheme 

The strongest constraints on mixing angles between ordinary fermions and mirror fermions arise 
from the conservation of e-, /i- and r- lepton numbers and from the absence of flavour changing 
neutral hadronic currents. In a particular scheme these constraints can be implemented at tree 



level [fL4|l . In this "monogamous mixing" scheme the structure of the mass matrix is such that 
there is a one-to-one correspondence between fermions and mirror fermions, due to the fact 
that the family structure of the mass matrix for mirror fermions is closely related to the one 
for ordinary fermions. 

Let us denote doublet indices by A = 1, 2, colour indices by c = 1, 2, 3 in such a way that the 
leptons belong to the fourth value of colour c = 4, and family indices by K = 1, 2, 3. Diagonal 
entries in the mass matrix for "normal" fermions will be denoted by an index ip, those for the 
mirror fermions by \ an d the off-diagonal mixing masses between fermions and mirror fermions 
by L or R (depending on the chirality of the "normal" fermion). In general the elements of the 
mass matrix for three mirror pairs of fermion families are diagonal in isospin and colour, hence 
they have the form 

f i (ip,xy,A 2 c 2 K2,AiCiK 1 — °A2A l C 2C l ^l 1 p jX y : K 2 K l ' 
^L;A 2 c 2 K 2 ,AiciK 1 = ^A 2 A 1 S C2 c 1 ^L^K 2 K 1 > ^R;A 2 c 2 K 2 ,A 1 c 1 K 1 = ^A 2 A 1 ^c 2 c 1 ^r-^kIk 1 ■ (1) 

The diagonalization of the mass matrix M can be achieved for given indices A and c by two 
6 ® 6 unitary matrices f£ and F R acting, respectively, on the L-handed and R-handed 
subspaces: 

F {Ac) \M^M) l F {Ac) , F R Ac) \M^M) R F R Ac) . (2) 

The main assumption of the "monogamous" mixing scheme is that in the family space 
fi^p, p, x , hl, fi R are hermitian and simultaneously diagonalizable, that is 

F (Ac) _ F (Ac) _ ( F^ \ 

Fl ~ Fr ~ I o F^ ) ' ( } 

where the block matrix acts in (ip, %)-space. The Kobayashi-Maskawa matrix of quarks is given 
by 

C = F i2c ^F^ , (4) 

independently for c = 1,2, 3. The corresponding matrix with c = 4 and A = 1 <-> 2 describes 
the mixing of neutrinos, if the Dirac-mass of the neutrinos is nonzero. (Majorana masses of the 
neutrinos are not considered here, but in principle, they can also be introduced.) 
An example for a mass matrix with "monogamous" mixing is the following: 

,(Ac) _UAc) n (Ac) , .(Ac) 

H'X'&Ki — A x r&K-iK! T~ °K 2 K 1 '-^ , 

(c) r r(c) (Ac) x (Ac) (Ac) <- x (Ac) / k n 

I^L;K 2 K\ - d K 2 K 1 L , VR;K 2 K x = X R VfaKiKi + d K 2 K 1 R , 

where ^ Ac \ A^ Ac \ 5^, A^ 1 ^, 5 R Ac ^ do not depend on the family index. 
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The full diagonalization of the mass matrix on the {ipL, ipR, Xl, Xr) basis of all three family 
pairs is achieved by the 96 ® 96 matrix 



/ (AcK) 

I cos a L 



n {LB) 



A'c'K'AcK 



Sa'aSc'cF^k ■ 







COS OL 



{AcK) 



(AcK) 

sin a L 







sin a 



(AcK) 



R 



sma 







cos a 



(AcK) 



(AcK) 
R 







V 







sin a 



(AcK) 







t R w COS Or 

M f M is diagonalized by (ZJJ)t M f MO iLR \ and MM* by C (i?L)t M M 1 {RL \ where {RL) is 
obtained from 0^ LR ^ by a^, <-> a R . 

In case of /xr = /il, which happens for instance in (|5|) if A# = and = 6l, the left-handed 
and right-handed mixing angles are the same: 



(AcK) 



(6) 



a (AcK) _ a JWO 



(7) 



In Ref. only this special case was considered. The importance of the left-right-asymmetric 



mixing was pointed out in Ref. ||17|| , where the constraints arising from the measured values of 
anomalous magnetic moments were determined. It turned out that for the electron and muon 
the upper limit is 

\a L a R \ < 0.0004 , (8) 
which is much stronger than the limits obtained from all other data Jl8| : 



a 



L) 



a 



I < 0.02 . 



* R ^ U.U- . (9) 

In case of the L-R asymmetric mixing the constraint (||) can be satisfied, for instance, if the 
right-handed mixing vanishes (or is very small): a R ~ 0. 

The hypothetical mirror fermions can be discovered at the next generation of high energy 
colliders. At HERA the first family mirror fermions can be produced via mixing to ordinary 
fermions up to masses of about 200 GeV, if the mixing angles are close to their present upper 
limits [[LJ3, BO]. At SSC and LHC mirror lepton pair production can be observed up to masses 



of a few hundred GeV [21]. This has the advantage of being essentially independent of the 
small mixing. At a high energy e + e~ collider, e. g. LEP-200 or NLC, mirror fermions can be 
pair produced and easily identified up to roughly half of the total energy, and also produced 
via mixing almost up to the total energy |22|] . 

In a model with fermion-mirror-fermion pairs the anomaly in the baryon- and lepton-number 
current is zero. Therefore in such models there is no anomalous fermion number violation at 
the elecroweak symmetry breaking scale. The baryon asymmetry of the Universe has to be 
produced at some higher energy scale, for instance, at the scale of grand unification [E3 . 



3 Fermion number anomaly on the lattice 



If mirror fermions do not exist, the lattice formulation of the Standard Model is problematic 
and the non-perturbative effects in the electroweak sector, like the anomalous fermion number 
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violation, cannot be studied by numerical simulations. There is, however, an approximation 
of the electroweak sector of the standard model which can be studied with standard lattice 
techniques, namely the limit when the SU(3) co i our (S>U(l)hypercharge gauge couplings are neglected 



3.1 Model for fermion number violation 

A simple prototype model is the standard SU(2)l Higgs model coupled to an even number 
2Nf of fermion doublets. In the standard model we have Nf = 6 (for simplicity, we consider 
Dirac-neutrinos, but the massless neutral right-handed neutrinos decouple [f24fl). One can take, 
for simplicity, Nf = 1 but the extension to Nf > 1 is trivial. The lattice action depends on the 
matrix scalar field (p x = 4>q x + i(p sx r s (with four real fields <fis=o,...,3) ari d the fermion doublet 
fields ^(1,2)*: 

S S sca i ar -\- Sf erm i 0n . (10) 
The standard scalar-gauge Higgs-model action is 



Sscalar = 4 E 4 m Tl (vlVx) + A [Tr {^ x f x )\ + ^ [ Tl (vlVx) ~ Tr (rf+A^M^)] > • (H) 



±1 

E 

The fermionic part contains the chiral gauge fields (with U XfX G SU(2) and Pl,r = (1 = F75)/2) 

U(L,R)xn = P(L,R)U xpL + P(R,L) , (12) 

and is given by 

S fermion = E { ^ [(^L eC ^lx) ~ {^lx eC ^2x) + (^2x eC Ai) ~ (l/j lx eC 'l/j 2x ) 



2 „ 



2x 



9 v bl^x^ix) ~ i^lx+^CU Lx ^ lx ) - (ipJ x eCip 2x ) + (^i x+fl eCU Lxa ip 2x ) 
+{^2x eC ^ix) ~ (^2x+fi u Rx^C^ x ) - ($ lx eCi? 2x ) + (^ lx+fl U Rx ^eC^ x ) 

+ (tPlRxGl^t^lLx) + (iplLxtpxGl^lKc) + ($%Rx G 2 ( Pt' t foLx) + ^ 2L x^xG 2 l\) 2 R x ) } • (13) 

Here e = ir 2 acts in isospin space, and C is the fermion charge conjugation matrix. The 
Yukawa- couplings G\ t2 can, in general, be arbitrary diagonal matrices in isospin space. In case 
of degenerate doublets G\ >2 are proportional to the unit matrix. 

Instead of the off-diagonal Majorana mass /xq and Majorana-like Wilson term (proportional 
to r), it is technically more convenient to consider a Dirac-like form with if) = ipi and the mirror 
fermion field \ defined by 

Xx = e^Ciplfr , Xx = ^2x eC ■ (14) 
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Figure 1: The value of the lattice sum on 100 4 lattice as a function af the bare mass fiQ. 
The upper curve is for periodic, the lower one for antiperiodic boundary conditions. 
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Figure 2: The same as fig. 1 on 200 4 lattice. 
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In terms of ip and \ one obtains the mirror fermion action for chiral gauge fields | I5| , which is 
well suited for studying the physically relevant phase with broken symmetry. 

In the symmetric (i. e. confinement) phase there is a natural alternative choice in terms of 
the reshuffled combinations 

i>Cx = ^Lx + XRx , ^Nx = XLx + i>Rx ■ (15) 

On this basis the vectorlike nature of the model becomes explicit (75's appear only in the 
Yukawa-couplings). The SU(2) gauge field couples only to the "charged field" ip C i an d the 
neutral doublet ipN has only its Yukawa-coupling. 

The fermion number is the difference of the number of fermions (?/>-field) and mirror fermions 
(x-field). The gauge invariant fermion number current can be defined as 



1 r 



Jxn — 2 



{^x+pTIhUlx^x) + (V^Tm^Iz/A+a) - {Xx+pTtvUitenXx) - (Xxl/JJlx^Xx+fi) ■ (16) 



3.2 Volume dependence 

The anomalous Ward-Takahashi identity can be derived, as usual, on a weak and smooth 
backgroud gauge field with field strength 

F^(x) = d,A s v (x) - cU* (*) + ge stu Al(x)A u u (x) . (17) 
In the continuum limit of the backgroud field the result for 2Nf fermion doublets is 

{dMx)), = lim^Ajj^}^- 4 = VW^WW.W) • ( 18 ) 
Here the lattice integral X is given by 



1 r n /ifc cos k\ cos k 2 cos k 3 cos fc 4 
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X{r, no) = j—^r / -=— — 3— [r^kl/ cos k a - fi k ]d 4 k , (19) 



(27T) 4 i- 7r (V+fi) 

and the notations are 

^ = 1^0+ 2 ^ 2 > K = sin k v > fc^ = 2 sin y . (20) 
The integral X is the same as the one occuring in the chiral anomaly, and one can prove (see 

e. g. mM) 

X(r, 0) = ^ 2 (independently from r) . (21) 

( |T5D and (pT|) show that the correct continuum anomaly is reproduced at vanishing bare 
(Majorana-) fermion mass fiQ = 0. It is, however, interesting to investigate the hq dependence 
of the lattice integral in ([19]). The numerical evaluation of the corresponding lattice sum X^ 
on finite (L 4 ) lattices shows that X = lini£_». 00 Z£ is very small, probably Z(/i ,r) = for every 
positive fiQ. (See fig. 1 and 2.) The deviation from zero for small /xq depends on the boundary 



condition and, as it can be seen, is a function of Lfj, . The value of Xl comes close to zero near 
Lfio ~ 10. This behaviour implies that at every positive /io the non-conservation of the current 
([16]) resulting from the topological charge of the background gauge field exactly cancels the 
non-conservation due to the finite bare mass. 

In case of a quantized gauge field, on a fluctuating background, the fermion bilinears and 
the topological density are renormalized and mixed with each other. The form of the anomaly 



equation becomes a matter of convention (see, for instance, |[271| ) 



4 Numerical simulations in a 2d U(l) model 

Before doing numerical simulations in the model for anomalous fermion number violation dis- 
cussed in the previous section, it is useful to study a corresponding U(l) toy model in two 
dimensions, which has often been studied in this context (see e. g. f28fl). 

The lattice action depending on the compact U(l) gauge field U XfM = exp(L4 M (x)), (/x = 1,2) 
and, for simplicity, fixed length Higgs scalar field <f)(x), \4>{x)\ = 1 can be written as 

S = /?E E [l-c 0S (F^(x))]-2«EE^(^ + A)^(^) , (22) 

where the lattice field strength is defined for // = 1, v = 2 as 

F^(x) = A u (x + p) - A u (x) - A^x + u) + A„(x) . (23) 
Real angular variables —tt < 6 Xfl < tt on the links can be introduced by 

U wli = exp^) , 6 Xfl = A„{x) - 2vr • NINT(A^x) /2n) . (24) 
Fermions in this two dimensional model are introduced in the mirror fermion basis (ip,x), 



according to (|T4T), but will not be explicitly considered here. 



4.1 Topological charge 

The topological charge of U(l) lattice gauge field configurations can be defined as a sum over 
the contributions of plaquettes. The basic assumption is the existence of a piecewise continuous 
interpolation of the gauge field |3D|. The gauge invariant topological charge on the torus 
corresponding to periodic boundary conditions is obtained either from the "transition functions" 
f29f or from the "sections" 0] of this interpolated gauge field. 
Introducing the plaquette angle — 7r < 9 Xflu < n by 

$x^iv = @xfj, @x+ji,v @xv @x+P,ft ^^^XflU i (25) 

one can show that the topological charge Q is given by 

Q = E «*ia • (26) 



S 



4.2 Gauge field kinks 



It follows from the definition in (|25f) that the integers n x \2 can have the values 0, ±1, ±2. The 
plaquettes with n x yi ^ can be imagined to carry "Dime strings" or "gauge field kinks" |3"2"|j . 
These local contributions to the topological charge play an important role in the dynamics of 
lattice U(l) gauge fields and has been intensively studied in the literature. For recent references 
in higher dimensions see, for instance, |33|, 



In the scalar-gauge model defined by the lattice action fl22|), where the pure gauge part 
corresponds to the Wilson action with compact U(l) variables, in a typical configuration there 
are a lot of gauge field kinks. As an example, a typical configuration at the large ft k values 
(3 = 8.0, k = 1.0 is shown in fig. 3. As one can see, a considerable part of the plaquettes has a 
kink (with positive n 12 ) or antikink (negative n 12 ). 

It is important to keep in mind that the individual terms in the sum over plaquettes ( |2T>| ) 
are not gauge invariant, in contrast to the total topological charge Q. Namely, performing large 
gauge transformations on the two ends of a link can create or annihilate a kink-antikink pair 
on the two plaquettes which both contain the link. 

In order to decrease the number of gauge kinks one can fix the gauge in some way. One 
possibility is to define the "minimal gauge" by minimizing the sum of squares (or absolute 



value) of the link angles |35 

E XX ■ (27) 

X fl=l 

This removes some part of the kinks in fig. 3, but still a lot of them remains. 

The role of the gauge kinks can be important for quantitites sensitive to the topological 
charge, or to the smoothness of the gauge field. Therefore it is instructive to study the model 
with some modified actions suppressing kinks. One possibility is to suppress the kinks by 
introducing a "chemical potential" for them p6| , which can in fact also completely remove 
them from the space of allowed configurations. 

Another modification of the gauge field action is suggested by the distribution of link angles 
in the minimal gauge. For small gauge couplings g 2 = <C 1, where most of the link 
variables are concentrated near 9 xtl = 0, there are, namely, also some secondary maxima at 
8 Xfl = ±7r/2. This is shown on the example of a configuration at (3 = 8.0, k = 0.8 in fig. 4. 
The secondary peaks are due to the fact that the Wilson U(l) gauge action is periodic for a 
shift of the plaquette angles by 4 ■ tt/2 = 2tc. Therefore, one way to push in the (3 — > oo limit 
all link angles to zero is to introduce a modified gauge field action like 

S 4 = 16ft E E I 1 - C0S K 9 ^ + 0*+A,»> - 0*» - ^W)/ 4 ]} • ( 28 ) 

x ^=l J i/=2 

The factor 16 in front of ft is introduced in order to have in the continuum limit of the action 
the same normalization as for ft 

The lattice action in (|28|) is not exactly gauge invariant because large gauge transformations 
of the links can cause a jump by 2n in the link angles, due to the limitation in (p4"|). This kind 



Figure 3: The plaquette contributions to the topological charge q = ni 2 on 64 2 lattice at 
= 8.0, k= 1.0. 
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link angles in minimal gauge 




of gauge non-invariance is, however, natural since we want to suppress exactly such large gauge 
transformations creating kink-antikink pairs. S4 remains gauge invariant if both link angles 
and gauge transformation angles are kept small, for instance, less than 7r/3 in absolute value. 
Therefore, the breaking of gauge invariance does not influence weak and smooth fields relevant 
in perturbation theory. 

The modified gauge action S4 suppressing kinks is by no means unique. One could, for 
instance, also use the "non- compact" formulation with the gauge field action 

Satc = §£ E F, u (x)F, u (x) . (29) 

We want, however, to stay close to the physically interesting model with SU(2) gauge field 
in four dimensions, where usually the compact formulation is used. The generalization of the 
gauge action S4 to SU(2) is possible and stays closer to the Wilson action than the non-compact 
formulation. 

The effect of reducing the density of gauge kinks by the modified gauge action 64 is dramatic. 
This can be seen on fig. 5. 
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Figure 5: The same as fig. 3 for the modified gauge action £4 in minimal gauge at /3 4 = 
2.0, k = 0.8. 
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topological susceptibility at beta=2. 
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Figure 6: The topological susceptibility as function of k at (3 = 2.0 on 64 2 lattice. 

4.3 Topological susceptibility and C-S numbers 

An interesting physical quantity is the topological susceptibility 

XQ = (Q 2 > - (Q> 2 . 



(30) 



This is strongly influenced by the density of gauge kinks, as is shown by figs. 6 and 7. Although 
the equality of bare couplings j3 and (3± does not mean that, for instance, the correlation lenghts 
are the same, it is clear that \Q gets substantially reduced by the suppression of gauge kinks. 
Even if (3 is increased to /3 = 8, xq is on ly reduced by roughly a factor of 10, not by a factor 
more than 1000 as between fig. 6 and 7. 

Another important topological feature is the distribution of Chern-Simons number. This 
can be defined on a L\ ■ L 2 lattice in the temporal gauge with U XuX2)2 = 1, x 2 = 0, 1, . . . , L 2 — 2 
by the sum of link angles 



I Li-l 

N C S(X 2 ) = — 0xi,X2,l 
2n x 1= 



(31) 



As an example, the distribution is shown in fig. 8 and 9. 

The time dependence (here x 2 dependence) of Ncs{x2) can be characterized by the expec- 
tation value of [Ncs{x 2 ) — N C s(x 2 — l)] 2 , which is illustrated by fig. 10. 

In summary, it seems that the density of gauge field kinks has a strong influence on the 
topological features in the two dimensional U(l) Higgs model. In order to clarify the physical 
significance of this, detailed studies has been started recently |J7 . 
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topological susceptibility at beta_4=2. 
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Figure 7: The same as fig. 6 with the modified action. 
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C-S numbers at beta=8. kappa=.8 
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Figure 9: The same as fig. 8 in larger bins. 
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Figure 10: The expectation value of the squared time difference of the Chern-Simons number 
on 64 2 lattice at /3 — 2.0 as a function of k. 
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